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> 1 INTRODUCTION 



ABSTRACT 

In the context of future space-based asteroseismic missions, we have studied the prob- 
lem of extracting the rotation speed and the rotation-axis inchnation of solar-hke 
stars from the expected data. We have focused on slow rotators (at most twice solar 
rotation speed), firstly because they constitute the most difficult case and secondly 
because some of the CoRoT main targets are expected to have slow rotation rates. 
Our study of the likelihood function has shown a correlation between the estimates of 
inclination of the rotation axis i and the rotational splitting Sv of the star. By using 
the parameters, i and 6i>* = Svsini, we propose and discuss new fitting strategies. 
Monte Carlo simulations have shown that we can extract a mean splitting and the 
rotation-axis inclination down to solar rotation rates. However, at the solar rotation 
rate we are not able to correctly recover the angle i although we are still able to 
measure a correct 61^* with a dispersion less than 40 nHz. 
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Understanding dynamical phenomena inside stars is one of 
the most important current challenges for stellar physics. 
During the last decades, helioseismology has allowed astro- 
physicists to constrain the internal structure and dynamics 
of the Sun. In the same way, asteroseismology will aim to 
improve our knowledge of stellar dynamics, especially con- 
vection and rotation. With future asteroseismic missions 
like CoRoT (Convection Rotation and planetary Transits, 
lBaglinll2003) . it will be possible for example to determine 
the extent of the convective region in stars and to extract 
information on rotation. Since 2003 the first Canadian satel- 
lite dedicated to asterosei smology, MOST (M icrovariability 
and Oscillations of STars. IWalker et alJl2003h . has been op- 
erational, beginning the space age for asteroseismology. 

Asteroseismology has already provided informa tion on 
the internal rotation of stars fe.g. lAerts et al.( |200j. for re- 
sults on a /3 Cepheid). However, the most accurate seismic 
information has been obtained for the Sun. Helioseismology 
has provided very accurate profiles of the internal rotation 
(see iThompson et al 1l2003|. and the references therein) a s 
deep as 0.2 Rg iCouvidat et alj|2003l : ICarcfa et ai]l2004al) , 
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thanks to the Solar Heliospheric Observatory and to ground- 
based networks. Because of the rotation of stars, modes are 
not single peaks but multiplets. The splitting of the mul- 
tiplet components gives information on the rotation speed 
in the acoustic cavity covered by the mode. Nowadays and 
in the near future, the asteroseismic observations will be 
limited to low-degree modes because of the absence of spa- 
tial resolution on the stellar surface. Thus, new inversion 
techniques have been developed and ch ecked to derive, for 
example, the radial rotation profi le (e.g. lGoupil et al]ll996l : 
'Loc hard. Samadi fc Goupif 2004) or to infer the latitudinal 
differential rotation (Gizon & Solanki 20J)4). 

Rotational splittings could be derived from the oscilla- 
tion spectrum along with the other mode parameters. How- 
ever, as we have learned from the solar case, the rotational 
splitting is harder to extract for low-degree modes because 
of the limited number of components in a multiplet. More- 
over, another difficulty appears in the stellar case: the angle 
of inclination (i) of the rotation axis, which determines the 
multiplet pattern, is generally unknown. 

iGizon fc Solankil (l2003h (hereafter GS03) have recently 
studied the simultaneous extraction of the splitting and the 
angle i from low-degree oscillation modes. We propose here 
to follow up their analysis by studying the potential of multi- 
mode fitting for more critical situations and by proposing 
automatable procedures. Our main objective is to determine 
the precision and limits in the determination of the rotation 
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of solar-like stars from a mission like CoRoT. To do so, we 
have simulated CoRoT-type observations (150-day long) of 
a Sun spinning at different speeds with different axis orien- 
tations. We have considered rather realistic signal-to-noise 
ratios (S/N) and we have focused on the particular situation 
of slow rotators (less than twice solar rotation). 

Our prelimin ary results have been outlined in 
iBallot et alJ i2004l) . The present paper fully develops this 
work. The layout of the rest of the paper is as follows. In 
Sect.|5|we describe the main properties of modes for a star 
under rotational effects. In Sect.|21we describe the techniques 
used to extract splittings and angle i from several modes 
together. In Sect.^Jwe present the results of our method ap- 
plied to several example cases. Finally, we discuss the fitting 
methods before concluding in the last section. 



2 OSCILLATION SPECTRUM OF A 
SPINNING STAR 

2.1 Mode properties 

Acoustic (p) modes in solar-like stars are excited by turbu- 
lent convective mo tions. Oscillations are damped but p erma- 
nently re-excited ijOoldreich. Murrav fc KximarlllOMl . The 
oscillation power spectrum of such modes can be modelled as 
a noisy Lorentzian profile. For a power spectrum classically 
computed with the Fourier transform of a regularly-sampled 
time series, this noise is a multiplicative exponential. A mode 
(n, £, m) - see below - is also characterized by its frequency, 
its amplitude and its FWHM. 

In solar-like stars, the width F of a p-mode depends only 
on its frequency v. For the Sun, the function V{u) shows a S- 
shape. There is a plateau in the range 2300-3200/iIIz around 
a value of 1 /iHz. At low frequency widths decrease rapidly 
and increase at high frequency (^e.g. iGarci'a et al. ,.2QQ4b.') . 

In the absence of rotation the frequency of a mode de- 
pends only on its radial order n and its degree £: we de- 
note it Uni- Modes are {21 -\- l)-times degenerate among the 
azimuthal order m. This degeneracy is removed by break- 
ing the spherical symmetry, especially by rotation. The fre- 
quency of mode (n, I, m) is expressed as Pnim ~ Vni + Svnim- 
The asymptotic first-order approximation, developed for a 
star spinning as a solid body with an angular veloci ty f2, 
gives 5vntm = —mSu with Su = Q/2tt I^Ledoux^ll951^ . We 
call Si/ rotational splitting (or simply splitting). 

For geometrical reasons, only low-degree modes have a 
sufficient amplitude to be visible in an oscillation spectrum 
due to the integration of the luminosity - or the radial veloc- 
ity - on the full stellar disk. Mode amplitudes also depend on 
their azimuthal order m. Calculations are rather straightfor- 
ward and can be found for example in GS03. Assuming the 
equipartition of energy between the different components of 
a multiplet {n,l), their amplitudes can be expressed as 

Anlm ~ aim{i)Vi anl = aim{i)Anl. (1) 

In this expression, the factor Vi is the mode visibility. It 
depends on the limb-darkening function, i.e. on the atmo- 
spheric properties. The visibility Vi decreases strongly when 
t increases: for (. — 1,..,5, we have calculated {Vt/Vo)^ = 
1.5, 0.53, 0.027, 0.0039, 0.00067, assuming an Eddington law 
for the limb-darkening function. For this reason, we expect 
to measure only modes ^ = 0, 1, 2 and probably a few (. = ?>. 
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Figure 1. An £ = 2 mode for three different speeds and angles. 

The factor apm(i) is the amplitude ratio of modes inside a 
multiplet. It is a purely geometrical term, depending on i, 
the angle between the line of sight and the rotation axis. 
This is true under only one condition, that the contribution 
of each stellar-surface element to the total fiux depends only 
on its distance to the disk centre. Even if it is not exactly 
true for velocity-fluctuation observations due to the rota- 
tion of the star (e.g. Hcimov 1999), this assumption stays 
very good for luminosity observations. The flnal factor of the 
mode amplitude ani ~ a{vni) depends mainly on the fre- 
quency and excitation mechanisms. We note A^i = Vla^i. 
This approach is valid for low rotation rates, when rotation 
can be interpreted as perturbation. 

Thus a mode (n, I) is modelled by a multiplet parame- 
trized by flve parameters (only three for ^ = 0): the central 
frequency u^t, the amplitude Ani, the width Tni common 
to all the components, the splitting 5v and the angle i. 

2.2 Classification depending on 

We have defined three different scenarios according to 5v: 

(1) 5v > Sinv , 

(2) V <5v< 5q2V , 

(3) Su<r ■ 

where 5o2V denotes the small separation i'n+i,e=o — i'n.e^2 
(around lO/xHz for the Sun in the range 2000-3000 pHz). In 
the first case, the components of different modes are mixed 
and it could be difficult to label each peak in a spectrum with 
the correct values of I, m and relative n. However, when this 
identification is done, all of the splittings Svnirn are accu- 
rately defined. In the second situation, mode identification 
does not pose any problem in general for good S/N and, as 
the components of a multiplet are well separated, splittings 
are easily measured. In the third and last case, the multi- 
plet components are blended. The effect on the amplitude 
ratio of a multiplet due to a given inclination axis is not 
always distinguishable from those of the splitting as illus- 
trated by Fig.0 For three different conflgurations chosen as 
an example, the mode profiles are nearly the same; only fine 
differences appear in the structure of profile tops. When an 
exponential multiplicative noise is taken into account, these 
differences are very difficult to catch. We have studied this 
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more challenging situation, corresponding to Su < l/iHz (for 
suns), i.e. Q. < 2^0 (r20/27r ^ 0.4 ^Hz). 



3 EXTRACTING THE MODE PARAMETERS 

3.1 Fitting modes: maximum likelihood 

Splittings and inclination angle should be deduced from 
the oscillation spectra at the same time as all the other 
mode parameters. For that, we use techniques developed 
and applied in full-disk-integrated helioseismology. Oscilla- 
tion spectra arc fitted with a maximum likelihood method as 
described by^Appourchaux, Gizon & Rabollo-Soarcs (1991). 
The power spectrum of a solar-like star is modelled as the 
sum of modes, modelled by multiplets, and a background 
noise, mainly due to convective motions (granulation, su- 
pergranulation) , and instrumental noises. The first step of 
the analysis is to remov e the backgrou nd, previously fitted 
following the model of iHarvevI 1^8^, to obtain a "flat" 
background. Then the modes are classically fitted alone or 
by pairs {£, £ + 2) according to the value of the small sepa- 
ration and the mode amplitudes. The residual background 
is considered as a constant inside the fitting window. As i is 
a new parameter compared to the classic helioseismic anal- 
ysis, we have explored its impact, especially on the splitting 
determination. 



3.2 Guessing and assumptions 

The fitting method needs guesses for the parameters to fit. 
This estimate is a starting point of the parameter-space ex- 
ploration by the algorithm maximising the likelihood. We 
denote by x the estimate of the parameter x. A crude esti- 
mation of the mode central frequency can be obtained, by 
looking for its centroid. The amplitudes and widths can be 
first determined on ^ = modes, which are insensitive to 
rotation. As amplitudes a„f and widths r„^ depend mainly 
on frequency, initial values for the modes 1^1 can be in- 
terpolated from those of £ = as follows: 



T _ ^2 T 7 _Vl An-1.0 + An,0 , s 

An-l,2 — J-^An.O, ^n-l,! — r—lT , (Z) 



rn-1,2 — r„_o, 



r„_i.o + r,i.o , , 

in-l,2 = T: • (■3) 



Determining the estimates i and Sv is not easy when 
multiplet components are not well separated. A first possi- 
bility is to fit each mode as a single Lorentzian. The compar- 
ison between the widths of two neighbouring modes £ — 2 
and £ = allows us to detect the presence of rotation (when 
j > 0), but a quantitative interpretation of this broadening 
is difficult because of the cumulative effects of i and Su. 

The sensitivity of the fitting to i and Sv has been tested 
along with the impact of the noise on i and 5u determination. 
Modes of interest have been fitted as follows: 

• Pairs £ = & 2 are fitted with eight parameters 
{Ao, A2, 1^0,1^2, r,b, 5i>,i): their amplitudes, their frequen- 
cies, a common width, the background level, the splitting 
and the inclination angle. Assuming r„_i,2 = F^.p is a good 
approximation as shown by the solar case JChapl in et alJ 




Figure 2. Likelihood function for one simulated spectrum in the 
plane (i, Sv) of the parameter space. All the other parameters are 
fixed to their simulated value. The power spectrum is taken in the 
range 2200-3000/iHz. The white colour corresponds to the highest 
likelihoods and the black to the lowest. The X is the simulated 
value (iojSuo) and the -|- is the maximum of the likelihood. The 
dashed line follows Susini = Suq sinio- 



'200(f). Thus, the parameter space is reduced as well as the 
computing time and the risk of non convergence. 

• Modes £ = 1 are fitted with six parameters 
{Ai, 1^1, Fi, 6, Sv, i). As the expected amplitudes of the ^ = 3 
modes are very small we do not fit them, although they are 
present in simulated spectra. Previous results have shown 
that such a simplification could introduce biases - especially 
to frequencies and splittings - if neglected modes are not 
sufficiently small and/or are too close to the fitted ones. We 
have been careful and we have verified that no significant 
bias has been introduced in our case. 

Different random values for l and 5u have been tested 
on several Monte Carlo realizations of the spectrum (see 
method in Sect. UHl . We have seen first that the solution 
is not unique and a certain dependence upon the first guess 
parameters is observed. For a given mode in a given real- 
ization, the fitting procedure can converge to some different 
couples (i, 5v) according to the initial values 1 and Sv. How- 
ever the main effect is due to the noise which has a strong 
impact on the estimation of (i, Sv) and disperses the results. 
Nevertheless, we observe a clear correlation in the determi- 
nation of both parameters. Results are organized along the 
curve: Svsini ~ constant = Svo sinio. We denote with an in- 
dex {Svo,sinio) the real (input) values of the parameters 
in the simulation. This can be explained by a study of the 
likelihood function for a simulated spectrum. Figure |21 shows 
such function in a plane (i, Sv) in the parameter space, with 
all other parameters fixed to their true value. We observe in 
such a plane a ridge following the curve Sv sin i = Svq sin io . 
Thus i and Sv are correlated. A new pair of independent 
parameters can be built: 



(i, Sv* 



with Sv* — Svsini. 



(4) 



Hereafter we use preferentially this new variable Sv* which 
is better suited than Sv for studying fitting issues and dis- 
cussing results. We did not find major differences between 
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using 5v* and 5v for the minimization routine, except for 
the error bars computed by Hessian-matrix inversion. 

3.3 Proposed strategy: multi-mode fitting 

With a classical fitting strategy, the determination of i seems 
very sensitive and tricky. We propose here another strategy 
aiming to improve the accuracy of the obtained value of 
i. We have fitted simultaneously several modes, as we can 
consider - to a first order approximation - that they have 
the same value of i and 5v* . 

(1) Choosing initial guesses. For A, V and v, see 
Sect. 13.21 We firstly fit pairs £ = & 2 and single modes 
1=1, using several (typically 20) different random values 
of t and 5u* . We use the average of all the obtained results 
as a better guess for these parameters. 

(2) Fitting simultaneously the modes (£=2,n — 1), 
(£=0, n) and {i=\,n). We used eleven parameters 
{A2, Aq, A\,V2,i'o,vi, ^0/2,^1, b,8u,i). We obtained also a 
series of values for {i,8u*). The mean values are noted (im 
and (m for mean value). We obtain in this way a first 
measurement of i and Q.. 

(3) Global fitting on a large range of the spectrum. 
Fitting several modes simultaneously, keeping free all the 
parameters, would be too costly in terms of computing time, 
and too delicate in terms of convergence. So we have decided 
to fix all the parameters but {i, 5u*) to their values deduced 
from the previous step. We choose as guesses i = im and 
Sv* = 5vm- Fitted results are denoted {ig,Sv*) {g as global). 



4 MONTE CARLO SIMULATIONS 
4.1 Defining the simulations 

The mode characteristics are derived from the observations 
of the Sun made by the GOLF instrument ( Global Oscilla- 
tions at Low Freauencv. iGabriel et al]ll99R^ . However the 
amplitudes have been adapted to simulate luminosity obser- 
vations instead of Doppler velocity measurements. We have 
treated the Sun as it was a main CoRoT target of magni- 
tude 6 observed during 150 days. For such a star, in the fre- 
quency range of interest (2200-3000 /^Hz), stellar noise have 
to dominate instrum ental and photon noise (see discussion 
inlMichel et alJl2005ll . Thus, S/N (as defined bv lLibbrech3 
Il992l) of the hightest component of a multiplet varies from 
15 to 150 for the ^ = 1 modes, from 4 to 45 for e = 2, and 
from 0.7 to 7 for £ = 3 (in configurations at 80°). The widths 
do not vary much (from 0.8 to 1.1 /iHz). 150-day (resolution 
ss 77 nHz) power spectra are created including ^ ^ 3 with a 
splitting Suq and an angle io that we want to simulate. In the 
chosen frequency range, there are six modes for each degree. 
This choice of interval results from a compromise: we have 
rejected modes with too low S/N (i.e. at low frequency) and 
peaks too broad, useless for our analysis (at higher frequen- 
cies). This will give us a lower limit of what we could obtain 
in the real case, with the hope that CoRoT will reach such 
modes. To introduce the no ise of each realization we follow 
iFierrv Fraillon et aP (ll998H by using a random exponential 
distribution which simulates the stochastic excitation. 

To test the analysis method a Monte Carlo simulation 



is done, i.e., we repeat A'^ times the method on the same the- 
oretical spectrum changing only the realization of noise. As 
the computing time required in each realization is quite high 
and we want to do many different cases, we have decided to 
limit the number of realizations A'^ to 100. The statistical sig- 
nificance of the results is small but it is enough to check the 
general trends of the solution. In order to verify our results 
we have increased A'^ to 1000 in some cases, e.g. io = 60°, 
Si/Q = 0.8 fiHz. The conclusions remain roughly the same. 
We have simulated six different configurations: two rotation 
rates fl = 1 and 2^0, i.e. Svq = 0.4 and 0.8/iHz, with three 
inclination angles io = 30, 60 and 80°. 

4.2 Star spinning twice as faster as the Sun 

This class of stars is the most favourable among those con- 
sidered. Results obtained with our strategy are satisfactory. 
A clear improvement is found, relative to classical fitting. 
The histograms of Fig.|21-a show the distributions of deduced 
parameters for each considered inclination angle. Both pa- 
rameter couples (imj^i'm) smd {igySv*) are plotted for every 
studied stellar orientation. We make three main comments: 

• in the three configurations, determinations of Sv* are 
non-biased and little spread: the dispersion is around 30nHz. 
Results given by averaging [Sv^) and by global fitting {5u*) 
are very similar. Global fitting does not lead to a noticeable 
change in this parameter in this situation. 

• On the other hand, the global fit (i.e. ig) brings, for 
i, a major improvement at low angle [i = 30°) according 
to averaged results im- Although there continue to be sev- 
eral highly spurious results (ig > 70°), a large number of 
realizations lie around 30°. 

• There is a slight bias on the i determination for the 
extreme values, but it remains smaller than the error bar. 

4.3 Star spinning as the Sun 

Fitting results for the configuration with Su = 0.4 fiRz are 
shown in Fig.|21-b. The study of the distributions of im, Sum, 
ig and SPg leads to two different conclusions for Si/* and i. 

• The Sv* distributions are quite narrow with disper- 
sions similar to the previous configurations (around 30- 
40 nHz). However a significant bias appears in the three 
cases, whereas it is negligible in the simulations at 2 JIq. 

• The angle i is not correctly extracted. The distributions 
are rather chaotic. However we have remarked that around 
a fifth of the realizations have given an angle of 90° . For the 
global fits of these low-splitting cases, this value behaves like 
an attractor during the likelihood-maximising process. 

We wanted to know if it is possible to extract the angle 
i from the selected modes in a configuration SI = SIq. To 
do so we have considered an idealized situation: we have 
performed "ideal" global fits. In such fits all the parameters 
- except i and Su* - are fixed to their exact values and not 
to the values deduced from a previous fitting step (cf. step 
#3 in the strategy 13.31 . Moreover the exact values io and 
Svq are chosen as guesses i and Su* . Thus all is optimized 
for fitting: only noise can influence the results. Results of 
this fitting method are plotted in Fig.j^b with dashed lines. 
Thus we can conclude that: 
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a) C = 2Cq (Suq = 0.8 ^tHz) 



b) C = (Suo = 0.4 ^tHz) 




0.3 0.4 
5v* (|iHz) 



Figure 3. a) Distribution of fitting results in the tliree configurations O = 2 Oq. On the left the angle i; on the right 5u* . Histograms 
plotted with solid lines show the results of global fits (ig and Su*); Histograms with dotted lines show the averaged results im and Su^. 
The dot-dash vertical line indicates the input value, b) Same as a) but for H = S^q. We have added distribution of the "idealized" global 
fits (of. text) plotted with dashed lines. 



• the bias on Su* disappears. It indicates that this bias 
was due to errors in the values to which the parameters were 
fixed. However, the dispersion stays the same: it is mainly 
generated by the noise. 

• the determination of i is not changed. Noise dominates 
above the signature of the angle and that seems inevitable 
in such data. 



5 DISCUSSION 

5.1 On the i/5!/-correlation and the law aim{i) 

In the framework of global helioseismology, IChapl in et al.l 
i200ll) have observed that changing amplitude ratios fixed 
inside the multiplets £ — 1,2 and 3 during the fit of solar 
spectra introduces a systematic bias on extracted splitting. 
We can understand the reason by studying the likelihood 
function shown in Fig. |21 Changing the amplitude ratio is 
similar to changing the angle i, thus it introduces a bias on 
splitting determination due to the correlation we have found. 
Our results generalise this observation. They show perfectly 
that we should be cautious of bias introduced when parame- 
ters are fixed, because of the c orrelation existing between the 
different parameters (see also lFierrv Fraillon et al]|l998ft . 

This analysis shows that it could be interesting to de- 
rive the angle i by other ways, like directly studying the light 
curve of stars an d trying to follow up modulations due to ac- 
tivity spots (e.g. iRucinski et al.l l2004'l . If such an additional 
constraint is available, the situation would become similar to 



the solar case and the amplitude ratios ae,n could be fixed a 
priori and individual splittings fitted. However the measure- 
ment of i must be sufficiently accurate (probably ~ 5-8°) 
otherwise the estimate of Sf will likely be biased. 

The results presented in this paper depend on the law 
we have used to link aim to i. For fitting, this law must 
be de fined a priori. As s hown once again by solar experi- 
ence JChaplin et alJl2004^ ■ when multiplet components are 
blended and not separated - which is the case here - fits are 
very sensitive to the chosen law atrn (i) • Luckily, for intensity 
observations, these ratios depend mainly on well-controlled 
geometrical considerations (cf. Sect. 12. ll . 



5.2 Limitation and improvements 

The situation can be improved if low-frequency modes are 
measured. For these modes the splitting can directly be mea- 
sured because of their finer widths. Then fixing the retrieved 
splitting can yield to a good estimation of i at higher fre- 
quencies where the multiplets are better defined and the 
infiuence of the stochastic excitation less important. 

In our simulations, we have assumed that the angle and 
the splitting are the same for all modes. While it is true 
that i is the same for every mode, Sv can vary for real stars, 
especially because of the differential rotation that could exist 
along the radius. However, for the Sun this variation is weak 
for low-degree modes in the studied frequency range. We 
could also attempt to extract not a mean splitting but a 
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results. In agreement with GS03, we find that at = 2 SIq 
we can retrieve both parameters in most of the cases, but 
with error bars improved by the global fitting, especially at 
low angle. However, at = JIq we have not been able to 
correctly recover the angle i. This result emphasizes the in- 
terest of having an independent measurement of the angle, 
but it has to be accurate enough to prevent the inclusion of 
a bias in the splitting determination. 
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Figure 4. Synthetic representation of biases and error bars for 
i and 6v* deduced from the simulations, in all the studied con- 
figurations. The crosses (x) mark the expected values (io,'5^'o). 
For f2 = 2 Cq cases, the boxes indicate the mean results and 
their dispersions. For f2 = 1 Q.q cases, only error bars on &v* are 
plotted because of the absence of good determinations of i. The 
two dashed lines are isorotations 5v = Suq = 0.4 and 0.8 /iHz. 

mean splitting for each degree, as was done for a first stage 
for the Sun fcf. lLazrek et al.lll99(ji . 

If £ = 3 modes have sufficiently high amplitudes in real 
observations to be correctly fitted, the results shown here 
would be improved. If they could be observed but with low 
S/N, we could tr y to use a so-called n-coUapsogram (cf. 
iBallot et al1l2004) to extract a mean splitting. This tech- 
nique can be summarized as follows: averaging the spectra 
of several i = 3 modes with different orders n, after remov- 
ing the £ = 1 neighbours, to enhance the S/N and define 
the multiplet better; and fitting the resulting spectrum. It 
needs a good determination of the central frequency for ev- 
ery mode, and small variations of Sv and the width (which 
is the case in the "plateau" frequency range). 

We can hope to derive even better results by denoising 
asteroseismic spectra. Filtering the spectrum and enhancing 
S/N cou ld improve the contra st of multiplets, guiding their 
analvsis. lLambert et alj i200 6h are proposing methods based 
on curvelet transforms permitting such denoising. 

This analysis performed on the Sun can be extended 
to sufficiently bright CoRoT targets with similar mass and 
evolutionary state. S/N will depend on the convective-noise 
level which will be observed in other stars. Some discussions 
on this topic ha ve taken place after the observations of Pro - 
cyon by MOST llMatthews et all2004l: bedding et all2005ll . 



6 CONCLUSION 

One of the challenges of present and future asteroseismic 
space missions is to extract stellar rotation rates and, wher- 
ever possible, the internal rotation profile. To do that, mode 
splittings Su have to be measured. We have studied the im- 
pact of the extra parameter i, appearing in asteroseismology, 
on the fitting. We have shown a correlation between Su and 
i, and defined a new parameter Sv* = Svsini. Strategies 
of multi-mode fitting have been developed, tested and vali- 
dated with Monte Carlo simulations. Figure 0] sums up the 
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